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ABSTRACT 

Context. We present the results from a combined study of the average X-ray spectral and timing properties of 14 nearby 
AGN. 

Aims. We investigate whether a "spectral-timing" AGN correlation exists, similar to the one observed in Cyg X- 
1, compare the two correlations, and constrain possible physical mechanisms responsible for the X-ray emission in 
compact, accreting objects. 

Methods. For 11 of the sources in the sample, we used all the available data from the RXTE archive, which were taken 
until the end of 2006. There are 7795 RXTE observations in total for these AGN, obtained over a period of ~ 7 — 11 
years. We extracted their 3-20 keV spectra and fitted them with a simple power-law model, modified by the presence of 
a Gaussian line (at 6.4 keV) and cold absorption, when necessary. We used the best-fit slopes to construct their sample 
distribution function, and we used the median of the distribution, and the mean of the best-fit slopes, which are above 
the 80th percentile of the distributions, to estimate the mean spectral slope of the objects. The latter estimate is more 
appropriate in the case when the energy spectra of the sources are significantly affected by absorption and/or refiection 
effects. We also used results from the literature to estimate the average spectral slope of the three remaining objects. 
Results. The AGN average spectral slopes are not correlated either with the black hole mass or the characteristic 
frequencies that were detected in the power spectra. They are positively correlated, though, with the characteristic 
frequency when normalised to the sources black hole mass. This correlation is similar to the spectral-timing correlation 
that has been observed in Cyg X-1, but not the same. 

Conclusions. The AGN spectral-timing correlation can be explained if we assume that the accretion rate determines both 
the average spectral slope and the characteristic time scales in these systems. The spectrum should steepen and the 
characteristic frequency should increase, proportionally, with increasing accretion rate. We also provide a quantitative 
expression between spectral slope and accretion rate. Thermal Comptonisation models are broadly consistent with our 
result, and can also explain the difference between the spectral-timing correlations in Cyg X-1 and AGN, but only if the 
ratio of the soft photons' luminosity to the power injected to the hot corona is proportionally related to the accretion 
rate. 
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1. Introduction 

The X-ray variability properties of AGN have been ex- 
tensively studied during the past twenty years. Significant 
progress has been achieved in the estimation of their X~ 
ray power spectral density functions (PSDs), which (among 
other things) can be helpful in the search for characteristic 
time scales in these objects. This progress has been made 
possible with the combined use of monitoring RXTE light 
curves (which are up to 5-10 years long in many cases) and 
shorter (1 to a few days long), high signal-to-noise, XMM- 
Newton and Chandra light curves. The results have shown 
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that the PSD has a ^ —2 power law shape at high frequen- 
cies and then, below a characteristic "break - frequency" , 
Vht, it flattens to a slope of ~ — 1 (e.g. Uttley, M'^Hardy, 
& Papadakis 2002; Markowitz et al. 2003, M'^Hardy et al. 
2004). Uttley & M'^Hardy (2005) (UM05 hereafter) list i^bf 
estimates for 14 nearby AGN. Using these results, M'^ Hardy 
et al. (2006) (M06 hereafter) demonstrated that the corre- 
sponding "break timescale", Tbr = l/j^bt, increases with 
increasing black hole mass, Mbh, and for a given Mbh, it 
decreases with increasing accretion rate, toe (in units of the 
Eddington limit). 

Knowledge of the X-ray properties of the Galactic black 
hole binaries (GBHs) has also advanced substantially dur- 
ing the past twenty years. Power spectral studies of Cyg X-1 
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in particular have been advanced considerably. Pottschmidt 
et al. (2003; hereafter P03) for example have used many 
RXTE observations between 1998 and 2001 to study the 
long-term evolution of the PSD. They also studied the en- 
ergy spectrum of the source and presented convincing ev- 
idence that its timing and spectral properties are closely 
linked: the characteristic time scales become shorter as 
the spectrum steepens. Axelsson, Borgonovo & Larsson 
(2006; hereafter A06), using several archival RXTE obser- 
vations, which cover all spectral states of the source, de- 
tected the same "spectral-timing properties" correlation as 
well. Shaposhnikov and Titarchuk (2006) also found that 
photon index and characteristic PSD frequencies are posi- 
tively correlated in Cyg X-1. 

The question whether AGN vary in a manner similar to 
that of GBHs is a long standing one. The availability of the 
better quality AGN light curves over the last few years (re- 
sulted from intense and extended RXTE monitoring cam- 
paigns to sample variability on very long time scales) has 
allowed a more quantitative comparison between AGN, Cyg 
X-1 and other GBHs (GRS1915-I-105 for example, in M06). 
In this work we used archival RXTE observations of 11 
AGN, together with data from the literature for 3 more 
objects, to estimate their average spectral slope and com- 
pare it with the characteristic frequencies that have been 
detected in their power spectra. We show that a positive 
"average spectral slope - characteristic frequencies" corre- 
lation exists, and we argue that this correlation is driven 
by accretion rate, for a given black hole (BH) mass: objects 
with a higher accretion rate relative to the Eddington limit 
should also have a steeper spectrum and a shorter charac- 
teristic time scale. We also compared the spectral-timing 
relation we found in AGN with a similar relation that has 
been detected in one GBH binary, namely Cyg X-1. The two 
relations differ, but by an amount that can be explained if 
we take into account the BH mass difference in AGN and 
Cyg X-1. Our results support the idea that AGN and GBHs 
vary in the same way. They also have interesting implica- 
tions regarding the nature of the X-ray source in AGN and 
GBHs. 



Table 1. Black hole mass and timing properties of the AGN 
in the sample. 



Name 


Mbh/IO" Mq 


Ref. 


/4,r/10"'Hz 


Ref. 


Fairall 9 


255±56 


P04 




MA03 


PG 0804-^761 


693±83 


P04 


9.6 


PA03 


NGC 3227 


42±21 


P04 




UM05 


NGC 3516 


43±15 


P04 


20+?« 


MA03 


NGC 3783 


30±5 


P04 




MA03 


NGC 4051 


1.9±0.8 


P04 


OUOU_3QQQ 


M04 


NGC 4151 


13±5 


P04 




MA03 


Mrk 766 


3.5 


WQ2 


6100+35°° 

DiUU_25oo 


V03 


NGC 4258 


39±1 


H99 


0.2+f 


MA05 


NGC 4395 


0.05±0.05 


V05 


19S00+''^°° 


V05 


MCG -6-30-15 


4.5 ± 1.5 


M05 


' "J-sio 


M05 


NGC 5506 


88 


PA04 


130+?f 


UM05 


NGC 5548 


67±3 


P04 


6.3ti« 


MA03 


Ark 564 


2.6±0.3 


B04 


230001^6600 


PA02 



Reference used for Mbh: B04: Botte et al. 2004; H99 - 
Bernstein et al. 1999; M05 - M'^Bardy et al. 2005; PAQ4 - 
Papadakis 2004; P04 - Peterson et al. 2004; V05: Vaughan et al. 
2005; W02 - Woo & Urry 2002. Reference used for i^br: MA03 
- Markowitz et al. 2003; M04 - M'^Bardy et al. 2004; PA02 - 
Papadakis et al. 2002; PA03 - Papadakis, Reig & Nandra 2003; 
V03 - Vaughan & Fabian 2003. 



Table 2. Summary of the i?^"!^? observations. 



Name 


Obs. Date 


No. of Obs. 


Fairall 9 


1996-11-03/2003-03-01 


672 


PG 0804-^761 


1999-01-24/2004-12-23 


259 


NGC 3227 


1996-11-19/2005-12-04 


1024 


NGC 3516 


1997-03-16/2006-10-13 


250 


NGC 3783 


1996-01-31/2006-11-08 


875 


NGC 4051 


1996-04-23/2006-10-01 


1268 


Mrk 766 


1997-03-05/2006-11-07 


221 


MCG -6-30-15 


1996-03-17/2006-12-24 


1216 


NGC 5506 


1996-03-17/2004-08-08 


627 


NGC 5548 


1996-05-05/2006-11-06 


866 


Ark 564 


1996-12-23/2003-03-04 


517 



2. Sample selection and data reduction. 

For the purposes of this study, we need to study AGNs with 
known v^f and average X-ray slopes. The sample of UM05 
is the best choice given the availability of the i/hf estimates, 
and the fact that these objects have been observed regu- 
larly with RXTE, over durations of at least a few years. 
This is more than ~ 1 — 2 orders of magnitude longer than 
the characteristic time scale of the sources in our sample. 
Most probably then, we have observed most of the possible 
flux and spectral variations that they exhibit. Furthermore, 
the study of the same RXTE observations that were used 
to estimate the power spectrum of the sources, offers us 
the possibility to determine their underlying spectral index 
around the same period of the power spectrum measure- 
ments. 

In Table 1 we list the 14 AGN in the UM05 sample, 
together with their Mbh and fbf estimates. For 11 of these 
objects we considered all the RXTE observations that were 
performed until the end of 2006; all data were taken from 
the public archive. Table 2 lists the date of the first and 
last observation, and the total number of RXTE observa- 
tions that we have used (second and third column, respec- 



tively). For the remaining three objects, namely NGC 4151, 
NGC 4258 and NGC 4395 we used data from literature 
to determine their average spectral shape, as explained in 
Sections 3.1, 3.2 and 3.3, below. 

We used data from the Proportional Counter Array 
(PGA; Jahoda et al. 1996) only. The typical duration of 
each observation was ~ 1 — 2 ksec. The data were reduced 
using FTOOLS v. 6. 3. The PGA data were screened according 
to the following criteria: the satellite was out of the South 
Atlantic Anomaly (SAA) for at least 30 min, the Earth ele- 
vation angle was > 10°, the offset from the nominal optical 
position was < 0°02, and the parameter ELECTRON-2 was 
< 0.1. Appropriate PC A background files Q, were used to 
calculate background model energy spectra in the 3-20 keV 
band. 



We used the latest background model for the faint ob- 
jects, pca_bkgd_cmf aintl7_eMv20051128 .mdl available from the 
RXTE Guest Observer Facility 
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3. Spectral analysis 

We extracted STANDARD-2 mode, 3-20 keV (where the 
PCA is most sensitive), Layer 1, energy spectra from PCU2 
only. After background subtraction, we used the XSPEC 
v.11.3.2 software package (Arnaud 1996) to fit a simple 
"power-law + Gaussian line" (to account for the Ka iron 
feature) to the spectrum of each observation. We used PCA 
response matrices and effective area curves created specif- 
ically for the individual observations by pcarsp. All spec- 
tra were rcbinncd using grppha so that each bin contained 
more than 15 photons for the statistic to be valid. 

The simple "power-law -|- Gausian" (hereafter PLG) 
model fitted well almost all of the individual RXTE spec- 
tra of the sources listed in Table 2. The only exception 
is NGC 5506, which hosts a heavily reddened Narrow 
Line Seyfert 1 nucleus (Nagar et al. 2002). Its central 
source is absorbed by neutral matter with column density 
of (3-4) X 10^2 cm-2 (Lamer, Uttley, & M'^Hardy 2000, 
Bianchi et al. 2003). For this reason we used a "wabsxPLG" 
model instead, with Nh fixed at 3 x 10^^ cm"^. This model 
fitted well most of the RXTE spectra of this source. 

3.1. NGC 4395 

NGC 4395 is a rather faint X-ray source, with an average 
2 — 10 keV flux of ^ 5 X 10"-^^ ergs cm"^ (Shih, Iwasawa 
& Fabian 2003), which has not been observed by RXTE. 
Shih et al. have estimated a mean spectral slope. Fobs, of 
1.46 ± 0.02 from the study of a ~ 7 day long ASCA obser- 
vation. Since Tbr ~ 500 s for this source (Vaughan et al. 
2005), wc can be reasonably confident that the ASCA ob- 
servation sampled most of the flux/spectral variations that 
the source displays. For this reason we adopted this value 
as the best estimate for its average spectral shape. 

3.2. NGC 4258 

The nucleus of NGC 4258 is heavily obscured in the X-ray 
band. Absorption column density estimates of the order of 
a fewxlO^'^ cm^2 have been reported in the past (Yang 
et al. 2007). The absorption also varies on time scales of 
months and years (Fruscione et al. 2005). NGC 4258 has 
been observed extensively by RXTE. There are 707 obser- 
vations (obtained until the end of 2006) in the archive. We 
reduced the data for all of them, and fitted the resulting 
spectra with a wabsxPLG model, keeping Nh fixed at 10^^ 
cm~2. The model fitted the RXTE spectra well. The mean 
spectral slope is equal to 2.02. 

Fruscione et al. (2005) found that the best-fit F values 
increase with decreasing instrumental angular resolution. 
They argued that this is due to the fact that a small ex- 
traction radius can isolate the central engine from the sur- 
rounding soft nuclear emission. It is unsurprising then that 
the mean spectral slope we estimated is similar to the best- 
fit value of ~ 2.1 that Fiore et al. (2001) reported from 
a study of BeppoSax data, and significantly steeper than 
the values obtained from the Chandra and XMM-Newton 
data analysis (Fruscione et al. 2005). For this reason, we 
used their results from the spectral analysis of 9 Chandra 
andXMM-Newton observations, which were performed be- 
tween May 8 2000 and May 22 2002. The weighted mean of 
their best fit spectral slopes is Fobs = 1.69, which we adopt 
as the average spectral slope estimate for this source. 



3.3. NGC 4151 

The X-ray spectrum of NGC 4151 is one of the most com- 
plex in AGNs, characterised by narrow and broad spectral 
features from soft to hard X-rays (e.g. Kraemer et al. 2005). 
The spectrum above 2 keV is afii'ected by absorption from 
a neutral absorber, which may even be patchy, and from 
warm material photionised by the central continuum (e.g. 
Beckmann ct al. 2005, de Rosa ct al. 2007). 

NGC 4151 has been observed regularly by RXTE. Data 
from 506 observations performed until the end of 2006 are 
stored in the archive. A simple PLG model yields Xicd ^ 
in 241 out of the 506 spectra. The use of a partial covering 
fraction absorption model, "pcf absxPLG", with the cov- 
ering fraction and Nh fixed at 0.55 and 1.8x10^^ cm"^ (de 
Rosa et al. 2007) improved the quality of the model fit but 
there are still more than a hundred spectra that the model 
could not fit well. Given the complexity of the NGC 4151 
spectrum and the low spectral resolution of PCA, it is not 
safe to use a more complex model to fit the RXTE data. 
Using the results of de Rosa et al. (2007) , based on the spec- 
tral study of 8 BeppoSax observations from January 1996 
to December 2001, we found that the weighted mean of the 
best fit spectral slopes is Fobs = 1-61. 



4. Results 

4.1. The mean spectral slope 

The top panel in Fig. 1 shows the sample distribution func- 
tion of the Mrk 766 best-fit spectral slope values. Fobs- 
A Gaussian fits well this distribution. This is the case 
with 5 more sources, namely: NGC 5548, PC 0804-1-761, 
NGC 3516, NGC 5506 and Ark 564. In the middle and bot- 
tom panels of the same figure we plot the Fobs distribution 
functions in the case of NGC 3227 and NGC 4051. An ex- 
tended tail towards low F values can be seen in both cases. 
In NGC 3227, these low-F values correspond to the tran- 
sient absorption by a gas cloud of column density in late 
2000 and early 2001 (Lamer, Uttley & M'^Hardy, 2003). As 
for NGC 4051, it is well established that it shows rather un- 
usual (among AGN) low flux states, which last from weeks 
to months, during which the X-ray spectrum becomes ex- 
tremely hard (Uttley et al. 2004, and references therein). 
Less pronounced tails towards low F values were also ob- 
served in the spectral slope distribution functions of Fairall 
9, NGC 3783 and MCG -06-30-15. 

Given the asymmetry in the spectral slope distributions 
we observe in some cases, we used the median of the_ dis- 
tribution as an estimate of the mean spectral slope. Fobs, 
for all sources in our sample. The results are listed in the 
second column of Table 3. The numbers in the parentheses 
correspond to the standard deviation of the Fobs distribu- 
tions. They indicate the scatter of the individual Fobs about 
their mean, and are representative of the typical range of 
the spectral slopes that we observe for each object. In the 
case of NGC 4258 and NGC 4151, the standard deviation is 
corrected for the contribution of the error of the individual 
best-fit r values, which in some cases was quite large. In 
the case of NGC 4395, the error listed is the one reported 
by Shih et al. 2003. 
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Fig. 1. The spectral slope distribution function of Mrk 766 
(top panel), NGC 3227 (middle panel) and NGC 4051 (bot- 
tom panel). 



Table 3. The mean spectral slope estimates. 



Name 


robs(o-r^^J 


r20%(crr2o%) 


Fairall 9 


1.82(0.16) 


2.01(0.08) 


PG 0804+761 


2.03(0.19) 


2.31(0.09) 


NGC 3227 


1.61(0.25) 


1.79(0.05) 


NGC 3516 


1.49(0.20) 


1.74(0.07) 


NGC 3783 


1.66(0.09) 


1.77(0.03) 


NGC 4051 


1.90(0.32) 


2.16(0.08) 


NGC 4151^ 


1.61(0.08) 


1.72(0.08) 


Mrk 766 


2.06(0.16) 


2.25(0.06) 


NGC 4258'' 


1.69(0.16) 


1.80(0.20) 


NGC 4395^^ 


1.46(0.05) 




MCG -6-30-15 


1.86(0.14) 


2.03(0.06) 


NGC 5506 


1.85(0.06) 


1.94(0.04) 


NGC 5548 


1.75(0.10) 


1.87(0.06) 


Ark 564 


2.51(0.22) 


2.82(0.16) 



Values in the parentheses correspond to the standard de- 
viation of the points about their mean. 

Values estimated using the results of de Rosa et al. (2007). 

Values estimated using the results of Fruscione et al. (2005). 
" Estimate taken from Shih et al. (2003). 



4.2. The "intrinsic" spectral shape of AGNs 

We found that most objects show significant spectral vari- 
ations (the results from a detailed analysis of these vari- 
ations will be presented elsewhere). If the Fobs variations 
correspond to variations of the intrinsic continuum slope, 
Pint, then Fobs will be representative of the average Fi„t 
as well. However, the AGN X-ray spectra, in some cases 
at least, are strongly affected by the presence of warm ab- 
sorbing material even at energies above 2 keV. For exam- 
ple, significant warm absorbing effects have been observed 
in NGC 3516 (Turner et al. 2005), NGC 3783 (Reeves et 
al. 2004), MCG -06-30-15 (Miller, Turner, & Reeves 2008) 
and Mrk 766 (Turner et al. 2007). In this case. Fobs will be 
a biased estimate of Fjnt . 



It is even possible that the spectral variations we ob- 
served are mainly caused by variations in a complex and 
multi-layered absorber, while Fjnt remains constant (e.g. 
Turner et al. 2007, Miller et al. 2008). If that is the case, 
then, for each individual spectrum. Fobs = Tint + AF, where 
AF < (as any absorption effects always result in flatter 
spectra). Since we expect AF to be different from one obser- 
vation to the other (due to changes in the covering factor, 
ionization state of the absorber etc), then fobs = Tint -I- AF, 
where AF is the mean of all the individual AFs, and should 
be negative, hence Fobs < Tint. 

It has also been suggested that the observed spectral 
variations are caused by the combination of a highly vari- 
able (in flux) power-law (with Fi„t=constant) and a con- 
stant reflection component (e.g., Taylor, Uttley & M'^Hardy 
2003, Ponti et al. 2006, Miniutti et al. 2007). Even in this 
case, we expect _that Fobs — Tint + AF (with AF < 0), hence 
Tobs = Tint -I- AF. 

The point is that, if the 2 — 20 keV spectrum is affected 
by absorption and/or reflection effects, then Fobs will be a 
biased estimator of Fint- One way to minimise the bias is 
to estimate the mean spectral slope using only the largest 
Tobs for each object since, in this case, AF will be minimum. 
For this reason, we estimated the mean of the Fobs which 
are larger than the 80th percentile of the distribution (i.e. 
the value below which 80 percent of the observations fall). 
In the case of NGC 4258, we used the mean of the three 
steepest F values of Fruscione et al. (2005) as an estimate of 
1^20% ■ Similarly, in the case of NGC 4151 we considered the 
mean of the two steepest F values of de Rosa et al. (2007). 

The 1^20% estimates are also listed in Table 3. The num- 
bers in the parentheses correspond to the standard devia- 
tion of the points in the 20% upper part of the Fobs dis- 
tributions, and indicate the scatter of these points about 
F2o%. The 80/20 dividing line is somehow arbitrary, and is 
mainly determined by the need to retain a sizable sample 
of Fobs values to estimate their mean in the case of sources 
with small number of observations. In any case, our results 
do not change significantly when we use the 90th or 70th 
percentile of the distribution for the sources with more than 
800 observations. Furthermore, the f 20% values are closer to 
spectral slope estimates which are resulted when complex, 
and more realistic, models are fitted to high quality X-ray 
spectra. For example, in Mrk 766, Miller et al. (2007) derive 
F = 2.38 ± 0.04 from the "principal component analysis" 
(fit to eigenvector 1). In NGC 3516, Turner et al. (2005) 
derive F = 1.82 ± 0.01 and 1.77 ± 0.02 from XMM-Newton 
observations at 2 different epochs. In MCG -6-30-15, the 
absorption model in Miller et al. (2008) derives a best-fit 
of F = 2.26 ± 0.02 (from the Suzaku data analysis) and 
F = 2.28 ± 0.01 from the study of the XMM-Newton long 
looks. In comparison, Minutti et al. (2007) in their blurred 
reflection model also derive a steep F = 2.26 ± 0.04 from 
the Suzaku data. Finally, in NGC 4051, in a recent paper by 
Terashima et al. (2008; PASJ, submitted), based on absorp- 
tion fits to Suzaku data, derive F = 2.04±0.01. The average 
difference between these spectral slope estimates and r2o% 
is less than '--^ 0.1. Therefore, 1^20% should be more repre- 
sentative of the intrinsic spectral slope of each source, if 
absorption and/or reflection effects are significant. 
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Table 4. The KendauU's t test results for the correlations 



plotted in Fig. 2. 


'''/PnuU 




r vs fbf 0.20/0.32 


0.21/0.30 


r vs BH mass 0.11/0.58 


0.10/0.62 


f vs l/norm 0.53/8.2 X 10"=* 


0.54/7.3 X 10"^ 



^Pnuii depends on both the size of the sample and the value of 



^ Results in the case we use Fon^ instead of F, 



4.3. The "spectral-timing" relation in AGN 

Filled squares in the upper panel of Fig. 2 indicate the 
"mean spectral slope vs characteristic frequency" relation 
for the AGN in our sample, using the fobs values listed 
in Table 3 and the i^bf estimates listed in Table 1. The 
plot in the middle panel shows the "mean spectral slope vs 
BH mass" relation, and in the bottom panel we plot fobs 
as a function of the characteristic frequency multiplied by 
the BH mass of each object (using the BH mass estimates 
listed in Table 1). We call the product j/bfxMBH as the 
"normalised characteristic frequency" , I'norm (in units of 
HzxMq). The crosses in the bottom panel of Fig. 2 indicate 
the "mean spectral slope vs t'norm" relation when we used 
^20% (instead of Fobs) as a measure of the mean spectral 
slope for the sources in our sample. 

Visual inspection of Fig. 2 suggests that, although the 
mean spectral slope does not correlate with BH mass, it 
may correlate positively with ly^f, and even more so with 
Vnorm' objects with steeper spectra have shorter character- 
istic frequencies as well. We reached the same conclusion 
even when we replaced fobs with f 20% in the first two panels 
of Fig. 2. To quantify the correlation of the mean spectral 
slope with i/'bf , BH mass, and i^norm, we used the Kendall's 
T test. The test was performed in the log-log space, and 
the results (r and P„uU, i-e. the probability that, under the 
hypothesis there is no correlation, r could be this large or 
larger just by chance) are listed in Table 4. We accepted 
a correlation to be "significant" if PnuU < 0.05. Values in 
third column of_Table 4 are the test results when we use 
f2o% instead of fobs- 

The Kendall's r results imply that only the "mean spec- 
tral slope vs Vnorm" Correlation is strong (i.e r > 0.5)_and 
significant. This is true irrespective of whether we use f 20% 
or fobs • The correlation remains significant even if we omit 
Ark 564 (i.e. the source with the softest spectrum and the 
highest characteristic frequency): r = 0.45, Pnuii = 0.03 for 
the f 2o%~i^norm relation^ and r = 0.46, Pnuii = 0.03 when 
we use Fobs instead of F2o%- We therefore conclude that, 
although the mean spectral slope does not correlate signifi- 
cantly either with power spectrum break frequency or with 
the BH mass, it does correlate positively with the break 
frequency when normalised to the BH mass of each object. 







^ We got the same results when we excluded from the sample 
the three sources for which the mean spectral slope estimation 
was not based on the analysis of a large number of RXTE ob- 
servations. For example we found that r = 0.38 in the case 
of the "r2o% - J^br" relation (Pnuii = 0.10), r = —0.13, and 
Pnuii = 0.58 in the case of the "f2o% ~ Mbh" relation, and 
r = 0.53, Pnuii = 0.02 in the case of the "^20% - i^norm" rela- 
tion. 



I 1 1 iiiii| — I 1 1 iiiii| — I 1 1 iiiii| — I 1 1 1 iiii| — I I 1 1 iiii| — I I 1 1 iiti 



4 



^^^^ 



I I I I I 

10' 

(1-17 1 

bi 

Tj 1 — I I I llll| 1 — I I I llll| 1 — I I I llll| 1 — I I I llll| 1 — I I I II 



10 10 10 10'^ 10 10 



+ 



I I I I I I I I LJ-LU 




1^ fl-nl^ I I I I I I I 



10 10 10" 10" 10 

^norJHz*M^) 

Fig. 2. Top panel: The "fobs - break frequency" relation 
for the AGN in our sample. Middle panel: The "Fobs - 
BH mass" relation for the same objects. Bottom panel: 
The AGN "mean spectral slope - normalised character- 
istic frequency" relation. Filled squares and crosses indi- 
cate the (Fobs, t'norm) and the (F2o%, ^'norm) data, respec- 
tively. The solid and dashed lines in the same panel indi- 
cate the best-fit power-law model to the (f obs, t'norm) and 
(r2o%:^'norm) data, respectively. 



To investigate whether a straight line or a power law 
model fits best the AGN spectral-timing relation, we con- 
sidered the 10 AGN that have well measured PSD breaks 
and we fitted their "average spectral slope - normalised fre- 
quency" data with a straight line in both the linear and log- 
log space, taking into account the errors on both the ^'norm 
and fobs (or F2o%) values. We found that the power-law 
model (i.e. a straight line in the log-log space) describes bet- 
ter the observed spectral-timing relation. For example, the 
linear model fit to the (f2o%j ^^norm) data resulted in a 
value of 25.1 for 8 degrees of freedom (dof). A linear model 
to the logarithm of the same data set resulted in 17.6 for 8 
dof. To quantify if this change in is significant, we com- 
puted the "ratio of likelihoods" , L1/L2 (Mushotzky 1982). 
It is defined as Li/L2=exp[(xi - X?)/2] =exp(AxV2), 
where Xi and X2 are the chisquare for the line fits in the 
logarithmic and linear spaces, respectively. We found that 
L1/L2 = 42.5. This result suggests that a power law is 
40 — 45 times more likely, than a straight line, to be the 
"correct" model for the spectral-timing data plotted in the 
bottom panel of Fig. 2. 
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To derive the best-fit power law parameter values, we 
used the "ordinary least squares bisector" method of Isobe 
et al. (1990) to fit the [log(fobs), log(fnoim)] data with 
a straight line of the form log(robs) = a + b\og{vnorm)- 
The best-fit parameter values are a — 0.09 ± 0.05 and 
b = 0.070 ± 0.015 (the best-fit model is indicated by the 
solid line in the bottom panel of Fig. 2). The dashed line 
in Fig. 2 indicates the best-fit straight line model to the 
[log(r2o%),log(^'norm)] data. The best-fit parameter values 
in this case are a = 0.11 ± 0.06 and b = 0.079 ± 0.019. 
In other words, a power-law model of the form Fobs — 
1.23(±0.14):yO„o^±o oi5 fits well the (f obs, I'norm) data, while 
in the case of the (r2o%, t'norm) data the best-fit power-law 
model is f2o% = 1.3(±0.2)fO„o^?±o.oi9^ ^^^^ ^^^^ 

though the best-fit slopes are small, they are significantly 
different than zero. The best-fit parameter values are con- 
sistent within the errors and their weighted mean value 
is a = 0.10 ± 0.04 and b = 0.073 ± 0.012. We there- 
fore conclude that the spectral-timing relation in AGN 
is well parametrised by a power-law model of the form: 
Fagn = 1.26(±0.12)z.So°S±ooi2. 

According to M06, the PSD break frequency depends 
on both the BH mass and accretion rate approximately 
as follows: i^hi oc toe/Mbh- Consequently, t'norm should 
depend on ttie only. Given this observational result, the 
Tagn oc i^ngm relation we found can be translated to a 
f AGN oc rii^-^' relation. In other words, our results suggest 
that the mean spectral slope in AGN correlates positively 
with riiE: objects with higher accretion rate should also have 
steeper spectra (on average). 

In the case of Comptonisation models, where thermal 
electrons in a corona above the disc upscatter soft photons 
emitted by the disc of temperature Tg, the produced X- 
rays have a power-law spectrum. When a soft photon of 
initial energy Cg is Compton scattered in the corona, it ac- 
quires an energy Aeg, on average, where A is the Compton 
amplification factor defined as A = (Ldiss + Ls)/Ls (ijiss 
and Ls are the power used to heat the corona and the inter- 
cepted soft luminosity, respectively). The produced X-rays 
have a power-law spectrum whose slope, F, depends on the 
temperature and optical depth of the corona. In general, 
one expects that when Lg increases, the corona cooling will 
be more efficient. Consequently, its temperature should de- 
crease and the resulting X-ray spectrum will be steeper (i.e. 
F will increase). 

Using the numerical Comptonisation code of Coppi 
(1999; eqpair in XSPEC), which is applicable in the case 
of a soft photon source located in the centre of a static, 
isothermal spherical corona, Beloborodov (1999) derives 
a relationship between F and Lg/Ldiss for different disc 
temperatures, Ts. For T, ^ few xW^K, which is appli- 
cable in AGN, he finds that Fagn ^ 2.?,{Ls/ Ldissf '^ ■ 
For low mass black holes like Cyg X-1, where Tg '^few 
xlQ^K, Fgbhs ~ 2.3(Lg/L<jjss)° ". In the case of a non- 
static, outflowing corona, Malzac, Beloborodov & Poutanen 
(2001) found similar results: Fagn ^ 2.2{Ls/ Ldissf'^'^ and 
Fgbhs ~ 2.2{Ls/ Ldiss)^'"^^ ■ 

We found that Fagn ^'iTor'n7 which, when combined 
with the t'norm OC riiE result of M06, implies that Fagn oc 
j^~o.o7^ Consequently, thermal Comptonisation model pre- 
dictions, i.e F cx {Ls/ Ldissf'''^'^ ™ Q-W ^ consistent with 
our results, i.e Fagn rhE° °^7 but only if Lg/Ldiss hiE. 



5. Comparison with Cyg X-1 

To compare the spectral-timing behaviour of AGN with 
that of Galactic black hole X-ray binary systems, we con- 
sidered the extensive spectral and timing observations of 
the best studied GBH, Cyg X-1. To maximise the spectral- 
timing range, we considered all available information, cov- 
ering "hard" state observations (from P03 and A06) and 
"soft" and "intermediate" states (from A06). 

It is important that a consistent measure of character- 
istic break-frequency is used for all the Cyg X-1 and AGN 
data. In the soft state the PSD of Cyg X-1 is well described 
by a power-law with only one "bend", with slope ~ — 1 at 
low frequencies and ~ — 2 at high frequencies. This shape is 
sometimes also parametrised as a power-law of slope -1 with 
an exponential cut-off (model 5 of A06). This shape also fits 
well the PSDs of almost all AGN that have been studied 
so far (except for Akn564, for which a double Lorentzian 
model flts best, M'^Hardy et al. 2007). For the soft state 
PSDs of Cyg X-1 therefore, we considered from model 5 of 
A06, the best-fit "turnover" or "bend" frequency at which 
the PSD slope bends from —1 to —2. 

In the hard state, the PSD of GBHs such as Cyg X-1 
can be fitted either as a doubly bending power-law of slope 

at the lowest frequencies, slope —1 at intermediate fre- 
quencies, and slope —2 at the highest frequencies or, where 
the signal-to-noise ratio is higher, as the sum of a number 
of Lorentzian-shaped components. Both P03 and A06 have 
opted for the second option, in modeling the Cyg X-1 PSD 
in the hard and intermediate state. To include therefore 
the hard state Cyg X-1 PSD observations, we must use the 
Lorentzian which, in the bending power-law parametriza- 
tion, is closest to the frequency at which the slopes change 
from —1 to —2. In the observations of P03, that Lorentzian 
is referred to as V2 , and has an average frequency around 2 
Hz. In the analysis of A06 (and also of Axelsson, Borgonovo 
and Larsson, 2005), the hard and intermediate state PSD is 
parametrised by the sum of two Lorentzians but also with 
the addition of a weak cut-off power-law. The Lorentzian 
with the higher frequency corresponds to V2 of P03 and so 
we used the frequency of that Lorentzian hcrc0 

In order to measure Fobs, P03 fit their spectra with a 
"power-law -I- multi-temperature disc-black body -I- reflec- 
tion" model from which the Fobs of the power-law can be 
taken. A06 fist (20-9)/ (4-2) keV hardness ratios (HR). We 
have used the empirical relationship between F and HR of 
Axelsson et al. (2005) to convert HR into Fobs- 

The resulting Cyg X-1 data are shown in Fig. 3. Open 
and filled circles indicate the (Fobs, ^'2) data for Cyg X- 

1 where V2 is the centroid frequency of the "second" 
Lorentzian in the P03 and A06 model fits to the Cyg 
X-1 PSD, respectively. These data define the Cyg X-1 
spectral-timing relation in its hard state. Crosses indicate 
the (Fobs, Vc) data for Cyg X-1, where Vc is the "bend" fre- 
quency of the "bending" power-law model which fits best 
the Cyg X-1 PSD in its soft state, according to A06 (their 
model 5). To convert the observed frequencies (both V2 and 

^ We note that if band-limited PSD power is fitted by 
Lorentzians, the Lorentzian at the upper band limit will lie at 
a slightly lower frequency than the bend frequency, if the same 
PSD is fitted by a bending power-law (cf Akn564, M'^Hardy et 
al. 2007). The difference, however, is typically only a factor of 

2 or less, much less than the range of frequencies covered here, 
and so we do not try and take account of it here. 
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Vc) to i^norm we assumed a black hole mass of 15 M©, in- 
termediate between the two published estimates of 10 Mq 
(Herrero et al. 1995) and 20 M© (Ziolkowski 2005). We can 
see that the Cyg X-1 hard and soft state data form a smooth 
continuous distribution in the spectral-timing plane shown 
in Fig. 3. 

Filled squares in the top and bottom panels of Fig. 3 
indicate the "spectral slope - characteristic frequency" re- 
lation for the AGN in our sample, when we use the Fobs and 
f 20% values, respectively. Although in both Cyg X-1 and 
AGN the characteristic frequencies increase as the energy 
spectrum steepens, Fig. 3 shows clearly that the respective 



"F 



relations are not the same. 



This difference can be explained by the fact that the ac- 
cretion disc is hotter in Cyg X-1 than in AGN. So, for exam- 
ple, in the case of an outflowing corona, as stated in Section 
4.3, F oc (L,/Ld,,,)0-07 in AGN but F oc {Ls/ Ld^ss)°■^^ in 
Cyg X-1. If (Ls/Ldiss) oc riiE, both in AGN and Cyg X-1, 
then we expect that Fagn oc rhg°^ and FcygX-i oc ihg^'^. 
Using the I'norm oc riiE relation of M06, these relations can 
be written as follows: Fagn oc v^'^J^ and FcygX-i oc t'norm- 



As a result, in the case when Fcygx 
expect that t-J 



Tagn, we should 



,0.07 

norm, AGN 



,0.55 
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^norm,CygX-l' l'^' ^norm,CygX-l 



^Qj^ (in the case of a static corona, we expect a similar 

relation, i.e. t'norm,CygX~i ^ ^^Sorm.AGN)- Open squares in 
both panels of Fig. 2 indicate the AGN data when we trans- 
form ^'norm.AGN to ^'nojfjj AGN' ^hc agreement now between 
the Cyg X-1 and AGN "spectral slope - characteristic fre- 
quency" relation is good (a more detailed derivation of the 
relation between t'noim,AGN and t'norm.CygX-i in the case 
of thermal Comptonisation models is presented in Section 
6.2). 



6. Discussion and conclusions 

We have used 7795 RXTE observations of 11 AGN, ob- 
tained over a period of ~ 7—11 years, to extract their 3-20 
keV spectra. We fitted them with a simple "power-law -I- 
Gaussian line" model, and we used the best-fit slopes to 
construct their sample distribution function. We used the 
median of the distributions, and the mean of the best-fit 
slopes which are above the 80th percentile of the distribu- 
tions, to estimate the mean spectral slope of the objects 
(the latter estimate is more appropriate in the case when 
the energy spectra of the sources are significantly affected 
by absorption and/or reflection effects). We also used re- 
sults from literature to estimate the average spectral slope 
of three more objects. The fourteen AGN that we consider 
in this work are nearby. X-ray bright objects, whose X- 
ray light curves have been studied in the past. Their PSDs 
have been accurately estimated, and characteristic "break 
frequencies" have been detected in them. 

When we combine the mean spectral slope estimates 
with the fbr estimates listed in UM05, wc find that: objects 
with steeper mean energy spectra have shorter characteris- 
tic time scales as well. This is the first time that such a 
"spectral - timing" correlation is detected in AGN. The re- 
sults we reported in Section 4.3 suggest that this spectral- 
timing relation in AGN can be parametrised as follows: 




10000 



Fig. 3. (Both panels) Filled and open circles indicate the 
spectral-timing data for Cyg X-1 in its hard state, using 
data from P03 and A06, respectively. Crosses indicate the 
spectral-timing data for Cyg X-1 in its soft state using data 
from A06. Filled squares indicate the "mean spectral slope 
- normalised characteristic frequency" data for the AGN in 
our sample, when we consider the Fobs and the r2o% val- 
ues listed in Table 3 (top and bottom panels, respectively). 
Open squares indicate the same AGN data when I'nornijAGN 



are shifted to 



,0.55 



AGN (^66 text for details). 



0.07 



6.1. The spectral-timing correlation in AGN 

The easiest way to explain this correlation is to assume 
that, for an AGN with a given BH mass, the accretion rate 
determines both the PSD characteristic frequencies (this 
has already been shown by MOB) and their energy spectral 
shape as well: the higher the accretion rate, the steeper 
the mean energy spectrum will be. Such a positive F — riiE 
relation in AGN has already been suggested/shown in the 
past (e.g. Porquet et al. 2004; Bian 2005; Shemmer et al. 
2006; Saez et al. 2008). A positive F — riiE relation has also 
been detected recently in 7 GBHs by Wu & Gu (2008), when 
they accrete at a rate higher than ^ 1% of the Eddington 
limit. Our results are consistent with the results reported 
in these papers. 

Using the above results and the M06 results, we cal- 
culate that fnorm(HzxMo) w SOOOriiE. When we replace 
i^norm in the relation of F « LSi^J^Q"^ we found in this work, 
we find that: F w 2.3mg°^. This formula between spec- 
tral slope and accretion rate could be helpful in deriving a 



8 



Papadakis et al.: AGN spectral and timing properties 



rough estimate for the accretion rate of the distant AGN 
that are detected in deep X-ray surveys, assuming that their 
variabiHty properties are similar to the properties of nearby 
AGN (this is supported by recent studies, see e.g. Papadakis 
et al. 2008). 

In the context of thermal Comptonisation models, riiE 
can affect the spectral slope F as it controls the strength 
of the soft disc photons, hence the cooling of the thermal 
plasma in the X-ray emitting corona. The greater the cool- 
ing by seed photons incident on the plasma, the softer the 
resulting X-ray power-law spectra are. Indeed, we found 
that Comptonisation models are consistent with the T — riiE 
relation that our results imply, but only if Lg/idiss cc inE. 

6.2. The comparison with Cyg X-1 

Both P03 and A06 have used data from RXTE observations 
that lasted for ~ 1 — 2 ks, and were performed every ~ 5—10 
days over a period of many years. If "time scales" scale with 
BH mass in accreting objects, then 1 — 2 ksec in Cyg X-1 
should correspond to a period of at least ~ 3 — 6 years in 
objects with Mbh > 10^ times the mass of the black hole in 
Cyg X-1 (like most AGN in our sample). This suggests that 
each one of the AGN points in Fig. 3 corresponds to just 
one of the Cyg X-1 points plotted in the same figure. We 
found that the AGN and Cyg X-1 "F — z/norm" relations are 
similar but not the same. In Section 5 we discussed briefly 
some implications of this result. In the paragraphs below 
we discuss the implications of our results in more detail. 

The main aim of the discussion below is to investigate 
the constrains that the AGN spectral-timing relation, and 
its comparison with the similar relation in Cyg X-1, impose 
on thermal Comptonisation models, based on the particu- 
lar assumption that both the spectral and timing prop- 
erties of accreting systems are driven by accretion rates 
variations. We point out that other interpretations are also 
possible; see for example Kylafis et al. (2008) for an al- 
ternative explanation of the Cyg X-1 spectral-timing rela- 
tion, which does not assume that X-rays are produced by 
thermal Comptonisation. We also point out that, given the 
small number of objects in our sample, and the unavoidable 
uncertainty in the derived parameters of the AGN spectral- 
timing relation, the values of the various parameters in the 
equations below are somehow uncertain. 

As we showed in Section 6.1, the AGN spectral-timing 
relation and the M06 results imply that Fagn ~ 2.3™^''^. If 
X-rays in AGN are produced by thermal Comptonisation, 
we expeclH that Fagn = 2.2(Ls/idiss)°'°^- Therefore, the 
observations [i.e. Fagn ~ 2.3TOg''^] are consistent with 
the thermal Comptonisation model predictions [Fagn = 
2.2(L3/Ldiss)°-°'], only if 

(Ls/idiss) « 2mE. (1) 

An obvious implication of this result is that, if a cer- 
tain fraction, say /s, of the total accretion power, Ptot 
(= TyTTiEC^, where 77 is the efficiency of the accretion pro- 
cess), is converted to disc luminosity (i.e. ig — /s^tot), 
while Ldiss = /diss-Ptot, then the ratio (/s//diss) should not 

* The discussion in these paragraphs is based on the assump- 
tion of an outflowing corona, hence we adopt the results of 
Malzac et al., 2001. Similar conclusions can also be drawn if 
we assume a static corona. 



remain constant for a given object, but should rather in- 
crease with increasing accretion rate. 

Suppose that X-rays from Cyg X-1 in its low/hard 
state (LH) are also produced by thermal Comptonisation. 
In this case, thermal Comptonisation models predict that 
Fcygx-i = 2.2(Ls/idiss)°■'^ or 

Fcygx-i « 2.4m°-i3, (2) 

if we accept that equation (1) holds in this case as well. 
According to Kording et al. (2007), the normalisation of the 
i^norm vs TOE relation for the hard-state GBHs is ~ 8 times 
smaller than the normalisation in the case of the AGN in 
our sample. So, if fnorm,CygX-i/LS ~ 375toe, as opposed to 

l^norm^AGN « 3000rhE, then TOE « 3 X 10~^J^norm.CygX-l/LS- 

We can now substitute toe in equation (2) to determine the 
F — fnorm relation for Cyg X-1 in LH state: 

TcygX-l/LS ~ (3) 

The dashed line in the top panel of Fig. 3 indicates this 
relation. The agreement between the predicted F — I'norm 
relation and the Cyg X-1 data is rather good up to i^norm ^ 
70. The discussion so far suggests the following picture: 

a) X-rays from the AGN we studied are produced 
by thermal Comptonisation. The F — fnorm relation we 
observe is consistent with the predictions of thermal 
Comptonisation models but only if the (Ls/-^diss) ratio, and 
hence the (/s//diss) ratio as well, increase proportionally 
with accretion rate. 

b) X-rays in Cyg X-1 are also produced by thermal 
Comptonisation. Taking into account the fact that the nor- 
malisation of the fc'norm — We relation is ^ 8 times smaller 
in Cyg X-1 than in the AGN in our sample (Kording et al. 
2007) , the predicted F — t'norm relation agrees well with the 
Cyg X-1 data up to fnorm ~ 70. 

c) In the case when Fagn = FcygX-i, we expect 

that (is/idiss)cygX-l ~ (^s/-C'diss)AGN' ^"^"^ "^AGN,E ~ 

"^CygX-i E (using equation 1). Consequently, toe.agn = 
™E CygX-i ^i^'i' since toe < 1, AGN should operate on 
a lower accretion rate than Cyg X-1 when the spectral 
slope is the same in both systems. Furthermore, since 

Fagn ~ l-S^^Som.AGN TcygX-l ~ l-l^^Sorm^CygX-l 

(when F < 2.1 — 2.2), then Fagn — Fcygx-i implies 

that 1.3t/^;„°^ AGN = l-l'^norni.CygX-l. ^ud I^norm.CygX- 1 « 

^•S'^norm AGN- Therefore, when the spectral slope is the 
same (and less than ~ 2.1 — 2.2) in AGN and Cyg X-1, the 
former should operate at a lower accretion rate but their 
characteristic time scales should be shorter than those in 
Cyg X-1 (when normalised to the respective BH mass), be- 
cause the normalisation of the AGN to — v-aom relation is 
significantly larger than the normalisation of the respective 
Cyg X-1 relation in LH state. 

The vertical dot-dashed line in the top panel of Fig. 3 
indicates the value z^norm 70. For Cyg X-1, this nor- 
malised frequency corresponds to t/i = 1 Hz and V2 ^ b 
Hz for the centroid frequency of the low and higher fre- 
quency Lorentzians, respectively. As fc'norm (i-c the accre- 
tion rate) increases even more, then i) the ratio v\lv2 in- 
creases as well (see Fig. 2 in A06), ii) the contribution of the 
Lorentzians to the root mean square variability amplitude 
decreases, and iii) the "bending" power-law component in 
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the PSD appears and its contribution to the source variabil- 
ity amplitude increases (see Fig. 7 in A06), i.e. the Cyg X-1 
power spectrum changes from a hard to a soft-state shape. 
At t'norm ~ 70, thc average observed spectral slope of Cyg 
X-1 is ^ 2.1 (see Fig. 3), while at higher frequencies the 
spectral slope is steeper. Consequently, the region defined 
by r < 2.1 and t'norm < 70 corresponds to thc "hard state 
region" for Cyg X-1 in the "spectral-timing plane" shown 
in Fig. 3. 

The discrepancy between the predicted spectral-timing 
relation and the Cyg X-1 data when F > 2.1 and i^norm ^ 70 
cannot be explained by the fact that the normalisation of 
the i^norm — "^E relation increases by a factor of ~ 8 in the 
high/soft state. If that were thc case, thc spectral-timing 
relation in this state should be similar to the one defined by 
equation (3), but with a sma/ier normalisation (opposite to 
what we observe). 

One possibility is that the (Lg/Ldiss) — relation 
(equation 1) changes in the high/soft state. However, in 
this case we would have to accept that the X-ray source 
does not operate in the same way in AGN and Cyg X-1: 
when i^norm > 70, both Cyg X-1 and thc AGN in our sam- 
ple follow the same Vnorm — relation (M06, Kording ct 
al. 2007). Therefore, as long as i^norm > 70, a given Unorm 
value implies the same accretion rate in both systems. The 
fact that the AGN spectral-timing relation is valid up to 
^'norm ~ 1000 should thcn imply that, for the same ac- 
cretion rate, the "(Ls/Ldiss)— accretion rate" relation is 
different in AGN and Cyg X-1. 

In Cyg X-1, r « 2.1 implies that 2.2(Ls/idiss)°-^^ « 2.1, 
and hence (Lg/idiss) ~ 0.7. Another possible explanation 
then for the discrepancy between the Cyg X-1 data and 
the predicted spectral-timing relation above fnorm = 70 is 
the following: equation (1) holds until (ig/i'diss) ~ 0.7, at 
which point the hot corona is significantly cooled down, 
and the thermal X-ray emission component is weak. It is 
possible then that at high accretion rates a separate, possi- 
bly non-thermal. X-ray component emerges, and dominates 
the X-ray emission in the soft state. If that is the case, the 
F — Ls/idiss relation we have assumed above is not valid, 
hence the predicted spectral-timing relation does not fit the 
data any more. 

Even if the picture drawn above is correct, there are 
important issues regarding the relation between AGN and 
GBH states which still remain unresolved. In particular, thc 
answer to the question whether the AGN in our sample are 
"soft" or "hard" state systems is far from clear. There are 
indications that they are analogous to Cyg X-1 in its soft 
state. For example, the radio emission in Cyg X-1 in its LH 
state is enhanced. On the other hand, Pancssa ct al. (2007) 
have shown that, for the same Lx/Le ratio, the radio lumi- 
nosity of Seyfert galaxies is ~ 8— 10 times lower than the ra- 
dio luminosity of hard-state GBHs, even when the BH mass 
difference is properly taken into account (7 of the objects in 
our sample are also included in their sample). Furthermore, 
the AGN in our sample follow the soft-state "characteristic 
time scale - accretion relation" in GBHs, and a "bending" 
power-law is the dominant component in the power spec- 
trum of those objects with good enough light curves to 
accurately study their PSD (e.g. NGC 4051, M'^Hardy et 
al. 2004; NGC 3227 and NGC 5506, UM05; MCG 6-30-15, 
M'^Hardy et al. 2005; and perhaps NGC 3783, Summons et 
al. 2007), implying a soft-state like PSD in these objects. 
High quality light curves for low accretion rate AGN are 



necessary to investigate whether any AGN with hard-state 
like power spectra exist or not. 

However, although the radio emission strength and the 
timing properties of many objects in our sample are soft- 
state like, their spectral properties are not, as the aver- 
age spectral slope is smaller than 2.1 in most cases. If in- 
deed the spectral hard-to-soft state transition corresponds 
to (Lg/Ldiss) ~ 0.7, at which point the thermal corona 
emission is weak, and a different component dominates thc 
X-ray emission, then we should expect this transition to 
happen when P > 2.15 in AGN. Given the AGN F - Vnorm 
relation, this slope corresponds to I'norm ~ 2500. At even 
higher normalised frequencies, we should then expect the 
AGN spectral-timing relation to break (like in Cyg X-1 at 
i^norm ~ 70). Obviously, more data are necessary to confirm 
that this is indeed the case in AGN. 

The data so far suggest that while in AGN the tim- 
ing properties transition from hard-to-soft state happens 
at least as low as Vnona ~ 100, the spectral properties tran- 
sition should happen at much higher accretion rates. This 
is opposite to what we observe in Cyg X-1, where both 
the spectral and timing properties change from the hard 
to soft state, at the same accretion rate (indicated by the 
value t'norm ~ 70). Perhaps the timing properties in ac- 
creting objects are determined by accretion disc variations 
only, and the hard-to-soft state transition materialises at a 
certain accretion rate, irrespective of whether the soft disc 
luminosity is strong enough for (Ls/-t'diss) > 0-7, i.e. strong 
enough to cool down the hot corona. In this case, we would 
expect the AGN hard-to-soft timing properties transition 
to appear at i/norm ~ 70 as well. Due to the cooler disc 
temperatm'c though, the AGN spectral soft-to-hard state 
transition happens at higher accretion rates (i.e. at a higher 
I'norm valuc in the spectral-timing plane of Fig. 3). Further 
progress in understanding the relation between AGN and 
GBHs can be made when we know how the accretion rate 
determines the characteristic frequency in accreting com- 
pact objects (assuming that it is just the accretion rate 
that determines I'norm in these systems). 
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